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Cold gas in the Intra Cluster Medium: implications for 
flow dynamics and powering optical nebulae 

Edward CD. Pope*, Thomas W. Hartquist, Julian M. Pittard 

School of Physics & Astronomy, University of Leeds, Leeds, UK, LS2 9JT 



00 

o 

O 

(N 

S3 
33 
>—> 

O 

Oh 
6 



> 

m 
m 

O 

oo 
o 



X 



20 June 2008 



ABSTRACT 

We show that the mechanical energy injection rate generated as the intra-cluster 
medium (ICM) flows around cold clouds may be sufficient to power the optical and near 
infra-red emission of nebulae observed in the central regions of a sample of seven galaxy 
clusters. The energy injection rate is extremely sensitive to the velocity difference 
between the ICM and cold clouds, which may help to explain why optical and infra- 
red luminosity is often larger than expected in systems containing AGNs. We also find 
that mass recycling is likely to be important for the dynamics of the ICM. This effect 
will be strongest in the central regions of clusters where there is more than enough 
cold gas for its evaporation to contribute significantly to the density of the hot phase. 

Key words: cooling flows, galaxies: active, galaxiesxlusters 



1 INTRODUCTION 

Line-emitting nebulae surround appr oximately a third of a ll 
Brightest Cluster Galaxies (BCGs) (jCrawford et alj|l999h . 
The occurance of these optical nebulae appears to have some 
correspondence to the short radiative cooling time of the 
cluster. Nearby examples, such as NGCI275, in Perseus, 
NGC4696 in Centaurus, and Af 795, show extended filamen- 
tary structures up to 50 kpc from the central galaxy (see 
lHatch et alj 120071 . and references therein) . Many of these 



BCGs also contain reservoirs of I0 8 - 10 11 ' 5 M of molecu- 
lar hydrogen (e.g. |Edg e| l200ll ). 

One of the prevalent emission lines in these nebu- 
lae is [NIIJA6584. Since [Nil] is collisionally excited by 
thermal electrons, this line measures the heating rate in 
the gas, whereas the Ha line measures the ionisation 
rate (jPonahue fc Void If99ll b A large heating rate com- 
pared to the ionisation rate leads to large [Nil] /Ha and 
[OII]/Hq ratios. To date, the power source remains un- 
known, although numerous heating processes have been in- 
vestigated. These include ionisation by the central active 
galactic nucleus (AGN), ionisation by hot young stars, heat- 
ing by X-rays from the ICM , conduction of heat from the 
ICM |Donahue et al.l l200Ch . and turbulent mixing layers 
jCrawford fc Fabianlll99a T Cosmic rays, preferentially dif- 
fusing along the magnetic field lines trailing behind rising 
bubbles, could possibly drive the excita tion in those fila- 
ment s that are located in bubble wakes jRuszkowsk i et all 
12007t h However, there are problems associated with each 
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possiblity listed above (see lHatch et al.ll2006h . As a conse- 
quence, the nature of the power source(s) remains unknown. 

Energy dissipation due to drag has, so far, not been 
considered as a potential source of energy for the filaments. 
In principle, dissipation by drag could be responsible for at 
least the heating of the optically emitting gas, and thus the 
[Nil] emission, if not the Ha emission. 

Recent numerical simulations suggest that the filamen- 
tary structures can be f ormed by an outfl ow stripping ma- 
terial from cold clouds l|Pope et al,ll2008h . These filaments 
only become long if the wind is sufficiently fast. The kine- 
matics and morphology also suggest that the filaments may 
be many tens of Myr old. Dissipation due to drag would 
necessarily occur all along the filament and would provide a 
spatially distributed source of energy that could power the 
optical emission. 

So far the presence of vast quantities of cold material 
has been ignored in most numerical simulations of the diffuse 
gas in galaxy clusters. This is partly due to the difficulty 
of implementing a cold phase in such simulations, but also 
because the importance of its presence has not been hitherto 
established. The effect of a cold phase may not be entirely 
negligible. In fact, the presence of significant quantities of 
cold gas might alter the flow dynamics in the central regions 
of galaxy clusters. This could occur through the addition of 
mass, stripped from cold clouds, or energy and momentum 
transfer associated with the mass loss from the clouds. 

In this article we estimate mass transfer parameters for 
a sample of galaxy clusters: Virgo, Perseus, Hydra, A2597, 
A2f 99, Af 795 and A478. We also investigate drag as a possi- 
ble power source for the optical emission of these nebulae. In 
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section 2 we discuss the constraints on cloud properties. In 
section 3 we concentrate on the significance of mass trans- 
fer between the phases. Section 4 focusses on mechanical 
energy as a power source for the optical emission, and we 
summarise in section 5. 



2 CONSTRAINTS ON COLD CLOUDS 

There is little direct evidence constraining cloud sizes in the 
ICM. For example, if the optical filaments are drawn from 
cold clouds the filaments are likely to be no wider than these 
clouds. However, the width of the filaments are usually con- 
strained to be < 1 kpc which is not a very strong limit. 
Maps presented bv lSalome et al.l (|2008T l show clumpy struc- 
tures in the CO emission. Given the spatial resolution of 
the technique, the upper limit on the sizes of these clumps 
is ~ 450 pc in the Perseus cluster. Such clumps may be com- 
plexes of giant molecular clouds, but it is difficult to tell. 

In the simplest cloud formation scenario there is a 
heat conduction length scale above which fluctuations 
are not stabilised and clo uds form, which is given by 
iBoehrineer fe Fabianl [jl989T l (following lFieldlll965h as 



A > 6 



0.01 



1/2 
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kpc, (1) 



where rj is the thermal conduction suppression factor, T is 
the ambient temperature, and n is the number density of 
the gas. Equation JTJ implies a minimum cloud mass of 
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(2) 

iMcKee fe Cowiel (|l977T ) present a similar argument. Equa- 
tion ((2]) indicates that smaller, lower mass clouds can form 
near the central regions of galaxy clusters, where the gas 
density is at its highest. More massive clouds may form fur- 
ther out, where the density is lower. 

However, the physics of the ICM is rather more com- 
plex than this description implies. Processes such as mag- 
netic turbulence are likely to play roles, with rj possibly 
as small as 10 -6 . Consequently, clouds much smaller in 
mass than 10 8 Mq may form. Given the difficulty of ob- 
serving single cold clouds, much of the work o n cloud prop- 
erties has relied on theoretical rea soning (e.g. iDaines et alJ 
1 19941 : lLoewenstein fe Fabian| [l990). However, there remains 
no clear picture of the mass and size of individual clouds. 
All we can say for certain is that, based on filament widths, 
and the CO data, the clouds are probably less than 0.5 kpc 
wide and there is likely to be a distribution of cloud sizes. 
Also, the cold mass should not be confused with the mass 
of optically-emitting ionised material, which can be deter- 
mined from the Ha emission and is much less than the total 
quantity of cold gas. 

If the most massive clouds are molecular, much of the 
cloud material may be at temperatures of tens of Kelvin. 
Pressure equilibrium with the surrounding X-ray emitting 
gas would imply a density of order 10 -20 — 10 -19 gem -3 . A 
spherical cloud of 10 s Mq would have a diameter of ~ 100 pc. 
Though the cloud mass is very uncertain, we will typically 
adopt values of 10~ 19 g cm -3 , 100 pc and 10 s M for the den- 
sity, diameter and mass where presenting numerical results. 



However, we will make the ways in which our results scale 
with the assumed values very clear. 

2.1 The motion of clouds 

The clouds will be referred to as the dispersed phase, while 
the intracluster gas will be referred to as the continuous 
phase denoted by subscripts d and c, respectively. 

Concentrating only on drag, we can write the equation 
of motion of a cold cloud through the ICM as 

dv 1 nD 2 
m ~dl = 2 D ^T Pc( u - V )\ u - V \> ( 3 ) 

where Co is the drag coefficient, D is the diameter of the 
cloud, p c is the density of the continuous phase, u is the 
speed of the continuous phase, and v is the speed of the 
dispersed phase. 

Dividing both sides of equation © by the cloud mass 
and assuming that the cloud is spherical with uniform den- 
sity, we find 

dv (it — v) 
dt t v 

where 



(4) 



4D 
3Co~ 



(5) 



Equations @ and © mean that the cloud's motion can be 
described in terms of the drag coefficient. The drag coeffi- 
cient is a function of the Reynolds number of the fluid flow- 
ing around an obstacle, and the Mach number (although for 
subsonic flows this effect is small). The Reynolds number is 
given by 



Re 



Dp c \u - 



(6) 



where fi c is the viscosity of the continuous phase. An ap- 
proximate relation for the subsonic drag coeffi cient for flow 
past a solid object at Re < 1000 is given by I Crowe et al.l 



24 4 
CD ~Re~ + R^7z' 



(7) 



So, Re = 1 gives Co = 28. At high Reynolds numbers 
(Re > 1000) the drag coefficient is approximately constant, 
Co » 0.45. This is not true if the Mach number is ^ 1 (e.g. 
lHendersonlll976l ). In such a circumstance the drag coeffi- 
cient increases. However, the effect of the Mach number is 
minimal for subsonic flows. 

Typical cloud velocities can be estimated by various 
means. In the Perseus cluster, the properties of the cold 
gas were probed with absorption of 21 cm ra dio emis- 
sion from a background radio source l|jaffeHl990h . Typical 
linewidths found were 100 — 500 km s -1 , which greatly ex- 
ceeds the expe cted therma l widt h. Similar linewidths were 
also found by lEded (|200lh and ISalome fe Combes! (|2003T ) 
from C O emission, and also f rom molecular hydrogen emis- 
sion bv lJaffe fc Bremer! (|l997h . These linewidths are compa- 
rable with inflow velocities in the ICM. However, this should 
not be taken as an indication that the clouds must be ap- 
proximately co-moving with the ICM. 

The Stokes number is a useful quantity to consider. It 
is defined as the ratio of the momentum transfer timescale 
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of the particle, t v , to the dynamic timescale of the fluid flow 
around the particle. The Stokes number is 



St 



Substitution for r v gives 



St = 



T v (u — v) 

D ' 



1 



(8) 



(9) 



which is independent of the cloud size and the relative veloc- 
ity. Since pd/pc 3> 1/Cb, we have St 3> 1. In this limit, the 
particle motion is not well coupled to that of the continuous 
phase and the two phases effectively move independently. 
This means that cold molecular cloud and the ICM veloci- 
ties are not coupled by viscosity. 

Subsonic bulk flows in the ICM are currently unde- 
tectable. Therefore, we need an alternative method to es- 
timate the quantity u — v that is im portant for our cal- 
culations. Thus, the IPope et al.l (|2008r i estimate of \u — v\ 
(derived from a model of the formation of filaments due to 
an AGN driven outflow) may provide a reasonable estimate 
for the Perseus cluster. 



2.2 Definitions 

Before calculating the mass transfer coefficients it is neces- 
sary to define some useful parameters. The first of these is 
the volume-averaged density of the dispersed phase (assum- 
ing identical clouds) 



Pd — ndm = Qdpd, 



(10) 



where rid is the number density of clouds of volume filling 
factor Qd and mass density pd- Similarly for the continuous 
phase we have 



p c — a c p c , 



(11) 



and the density of the mixture can be written p m = p c + Pd- 
However, for ease of comparing with observations, it 
makes sense to write equation (|10[) in a different way. The 
volume-averaged density is simply the total mass of the dis- 
persed phase divided by the total volume of the system 



Pd = 



Md 

v • 



(12) 



where V is the volume of a spherical shell, and Md is the cold 
mass within the shell. Similarly for the continuous phase we 
have 



Pc = 



3 MASS EXCHANGE 



Mc 
V 



(13) 



It is possible to assess the extent of the coupling between 
the phases by eva luating so-called 'coupling' parameters 
|Crowe et al.lll998l ). The continuous phase is described by 
the density, the temperature and the velocity field. The par- 
ticle/dispersed phase is described by the concentration of 
particles, their size, their temperature, and their velocity 
field. Coupling can take place through mass, momentum and 
energy transfer between phases. Mass coupling is the addi- 
tion of mass through evaporation, or the removal of mass 



from the continuous phase through condensation. Momen- 
tum coupling is the result of the drag force on the dispersed 
and continuous phases. It can also occur due to mass trans- 
fer between the phases. Energy coupling occurs through heat 
transfer between the phases as well as the dissipation of ki- 
netic energy due to drag. Thermal and kinetic energy can 
also be transferred between the phases by mass transfer. A 
quick numerical exercise demonstrates that estimates of the 
energy and momentum coupling, except that associated with 
mass exchange, are generally subject to large uncertainties, 
and therefore unreliable. 



3.1 Definition of the mass coupling parameter 

Suppose that there are rid clouds per unit volume in a cubic 
box with side X. If each cloud evaporates at a rate m, the 
rate at which mass is injected by the dispersed phase in this 
volume is 



M d = n d X rh = Nrh. 



(14) 



The mass flux of the continuous phase through this volume 
is 



M c ~ p c uX 2 , 



(15) 



and t he mass coupling parameter is defined as (jCrowe et al.l 
1998) 



Md 
Mc 



(16) 



If II mass -C 1, then the effect of mass addition to the con- 
tinuous phase would be insignificant. If II mass 3> 1, then 
the mass added to the system would dominate. The latter 
would be unphysical in galaxy clusters if it occurred globally, 
since it would imply that the total mass of gas was not con- 
served. However, in a restricted region this scenario would 
be consistent with a fountain flow, where the clouds are re- 
plenished over and over again. In the case where n mass ~ 1 
both condensation and evaporation are important. 

The mass coupling parameter is likely to be a function of 
radius in galaxy clusters. Therefore, equation (|16p is of little 
practical use, since we cannot compare like with like: the 
mass flow rate of the continuous phase must be evaluated at 
a single radius, while the mass injected into the continuous 
phase corresponds to a volume. 

In galaxy clusters, Md(r) is the net rate of mass added 
to the flow from the clouds (evaporation - condensation), 
while M c (r) is the mass flux of the hot gas. Several scenar- 
ios are possible. If there is no evaporation/condensation (i.e. 
no clouds at all), then M c = M x , where M^ is the pure 
flow rate of hot material. The accretion rate of new mate- 
rial onto the central galaxy, M ga i, is given by the value of 
the pure flow rate at the characteristic radius of the galaxy: 
Moo(rg a i). If there is condensation of material out of the hot 
phase, but no subsequent evaporation, then M c < M^. Of 
course, if the clouds also accrete onto the central galaxy, 
M ga i may still be equal to M^. A third possibility is that 
material evaporated from the clouds may contribute a sig- 
nificant fraction of M c , and so M c may exceed the actual ac- 
cretion rate of "new" material from infinity (i.e. M c > Moo). 
Such a scenario would require that the mass in the clouds is 
predominantly from material that was previously accreted 



© 2008 RAS, MNRAS 000, [jM] 



4 E.C.D. Pope, T.W. Hartqmst, J.M. Pittard 



onto the central galaxy, and which has at some later time 
been thrown out of the central galaxy back into the ICM, 
perhaps through AGN activity. 

Rates of mass deposition from the continuous phase 
to the cloud phase are inferred from X-ray data. The 
mass deposition in a typical cluster is consistent with a 
mass flow rate that increases li nearly with radiu s, Mobs = 
-Mo(r/r coo i), at r < r cool (e.g. iPope et ~all l2006l ). Here M 
is a constant and r coo i is the radius at which the ICM's radia- 
tive cooling time equals the Hubble time. If no evaporation 
of clouds occurred it would be reasonable to assume that the 
flow rate of hot gas at r coo i is M c — Mo and that the flow 
rate of hot gas at r < r coo \ is M c = M$r /r coo \. Now imagine 
the scenario where for r d < r < r coo i there is pure conden- 
sation, while for r < r d material may also evaporate from 
clouds back into the hot phase so that there is a complex 
exchange of mass between the hot gas and cold clouds (i.e. 
condensation and evaporation, and re-condensation and re- 
evaporation) . Much of the condensing gas at r < r d may be 
gas that previously evaporated from clouds recently expelled 
from the central galaxy in a fountain flow, or via AGN activ- 
ity, for example, rather than material that fell from r ^> r^- 
Thus, the instantaneous accretion rate onto the galaxy may 
exceed the long-term average (i.e. M ga i M 00 (r ga i)). There- 
fore, the observationally inferred mass flow rate Mo may ac- 
tually be an upper bound on the pure mass flow rate, Moo, 
with the real flow rate of "new" material being [3 Mo. The 
flow rate of new material in the hot phase is 



M = /3(r)M 



f cool 



(17) 



At r > r d , (3(r) = 1, while for r < r d , < (3(r) < 1. Let 
/? be the average value of /3 integrated between < r < r d , 
If there is not much recycling of material, then /3 will be 
close to unity. The average mass flux of new material in the 
hot phase at r < r d is then f3Mor d /r coo \. Using this, and by 
setting m/m = r m as the average ablation timescale of any 
cloud, we can write 



1 



^"cool 



( M d 



\/3Mo 



(18) 



If (3 is close to unity, then n mass is small everywhere. Al- 
ternatively, if P is significantly less than unity, n mass could 
be close to unity or larger. Equation (|18p is written in this 
way because we wish to express it in terms of observable 
quantities, such as the total mass of cold clouds, Md. r c \ is 
likely to be approximately equal to the maximum distance 
of clouds from the central galaxy. This scenario is pictured 
schematically in figure 1. 

3.2 Estimated mass-coupling parameter values 

We conside r the s ample of seven galaxy clusters studied by 
IPope et all (|2006h . Mo, r CO oi and Md are the mass flow rate, 
the cooling radius and the mass of cloud gas, respectively. 
Values for them and for the mass coupling parameter (based 
on the following discussion) are given in table 1. 

The rate at which mass is ablated from a cloud, m, 
might be estimated by the rate at which momentum associ- 
ated with the wind is transferred to the cloud, divided by a 
speed, c Sja . This speed would take a value anywhere between 
the sound speed of material in the cloud, and the difference 




Figure 1. Schematic view of the mass cycle described in the 
text. Material flows inwards through r coo i and a fraction might 
condense, so, for r d < r < r coo i we have M c (r) < Moo. For r < r d 
condensation and evaporation processes both occur as well as re- 
condensation and re-evaporation. 



Table 1 . Parameters for the sample of galaxy clusters used in 
deriving the cold gas masses and coupling parameters. Values 
in columns 2-3 are taken from Pope et al. (2006) and references 
therein. The values in column 4 are taken from Simionescu et 
al. (2007) for Virgo and from Edge (2001) for all others. Col- 
umn 5 shows the mass coupling parameter for a region enclos- 
ing 30 kpc around the cluster centre. The values in column 5 
could be greater if /3 < 1. 



Object 


M (M o yr" 1 ) 


f'cool (kpc) 


M d (1O 9 M ) 


n m ass 


Virgo 


1.8 


35 


0.5 


l 


Perseus 


54 


102 


17 


0.1 


Hydra 


14 


100 


<2 


<0.4 


A2597 


59 


129 


8.1 


0.1 


A1795 


18 


137 


<2.7 


<0.4 


A2199 


2 


113 


<2.9 


<3 


A478 


150 


150 


4.5 


0.06 



between the cloud speed and the average speed of material 
in the filament drawn from the cloud. The appropriate value 
of c Sja depends on how effective the viscous coupling (due 
to turbulence or any other mechanism) between the cloud 
material and the wind material is. The rate at which m ass 
is lost by the cloud would then be (e.g. IPope et al.ll2008l ) 

_ vi \2 

m = Ap c 



(u 



(19) 



where A is the cloud cross-sectional area. The ablation 
timescale is 



m 



Tm = 



'<•() 



(20) 



m \ p c J \u — v J \u — v t 
where ro is the initial cloud radius. We are aware of the 
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Table 2. Cold gas masses within 30 kpc calculated by setting 
n maS s = 1, ft = 1 and assuming r m = 10 9 yrs. The observed 
mass for Virgo is taken from Simionescu et al.(2007), all others 
are taken from Edge (2001). 



Object 


A/ d /10 9 M 


Observed mass/10 9 Mq 


Virgo 


1.5 


0.5 


Perseus 


15.9 


17 


Hydra 


4.2 


<2 


A2597 


13.7 


8.1 


A1795 


3.9 


<2.7 


A2199 


0.53 


<2.9 


A478 


30 


4.5 



large uncertainties in the parameters, but consider the fol- 
lowing estimates suggestive: ro ~ 50 pc, p&j p c ~ 10 6 , 
u — v ~ 10 8 cms _1 , c Sia = 10 6 cms _1 and r m ~ 10 9 yrs. 
Finally, we take the outer radius of cold material to be 
rd ~ 30 kpc. The mass coupling parameters calculated from 
equation (|18|) are shown in column 5 of table 1. These values 
suggest that in the central regions at least, the mass coupling 
may be a two-way process in some clusters, and mass-loss 
from clouds may be significant in comparison with the mass 
deposited within the central regions due to the actual flow 
of new material. Note that the values in table 1 are based 
on the highly conservative assumption that [3 = 1. If f3 < 1, 
then Il mass could be significantly higher. A caveat is that the 
estimated value of n maaa also depends on the actual mass 
loss rates from the clouds, which are very uncertain. Nev- 
ertheless, it seems that the mass coupling is probably an 
important consideration that has been omitted in numerical 
simulations of the ICM to date. 

The mass coupling parameter may also be interpreted 
in a slightly different way. The inflow of material, due to 
cooling, will deposit material onto the central galaxy of the 
cluster. Current thinking suggests that AGN activity lifts 
a fraction of this back into the ICM - this allows II mass to 
exceed unity. We can estimate the mass of cold gas in each 
cluster, by setting n mass = 1, and rearranging for Md = 
/3r m (30kpc/r coo i)Mo- These values are given in table 2, for 
{3=1. The estimated mass of cold material will be less if the 
ablation timescale of the clouds is smaller than assumed. 

The values in column 2 of table 2 compare relatively well 
with the observed quantities of cold gas in column 3. This 
suggests that the ablation timescale given by equation (|20p 
may be a reasonably accurate description of the process. 
Furthermore, since r m oc (it — v)~ 2 , it also suggests that 
it — v cannot be too different, on average, to 10 s cms -1 . The 
small differences between the values may be caused by slight 
differences between the assumed and actual value of r m , by 
recent AGN outbursts which may have temporarily lifted 
significant quantities of cold material, or by differences in 
the cloud masses. 



4 POWERING THE OPTICAL EMISSION IN 
GALAXY CLUSTERS 

In clusters of galaxies, the momentum transfer and dissi- 
pation of energy by drag into the hot gas are generally 
insignificant compared to gravity and radiative cooling in 



the hot gas. However, energy is still dissipated by drag in 
significant quantities. This is interesting because the ra- 
tios of H2 ne ar and mid-infrared to H recombination line 
strengths (e.g. lJaffe et"aT]|200ll . l200l I Johnstone et al.ll2007r ) 
are co mpatible with heati ng being due to dissipative pro- 
cesses ()Ferland et al.l l2008h . It is possible that energy dissi- 
pation due to drag constitutes such a process. 

We have already described ho w CO maps show clumpy 
structure (e.g. ISalome et al.| [2008) and that the upper limit 
on the si zes of these c l umps , in the Perseus cluster, is 
~ 450 pc. Sal ome et al.l |2008h suggest it is possible that 
these clumps are complexes of giant molecular clouds. If 
this is true, the limitations of spatial resolution mean it is 
impossible to trace the emission from a single cloud. Conse- 
quently, theoretical estimates of the emission due to physi- 
cal processes involving the interaction between hot and cold 
phases, must be applied to a volume containing multiple 
clouds, rather than to an individual cloud. The total me- 
chanical energy injection rate required to produce filaments 
of average length I, with a combined mass Md and a velocity 
shear along their length of Av, is 

Av 3 

E = SM d ^j—, (21) 

where I is the typical length of the filaments and 8 is a factor 
that takes into account the distribution of mass in terms 
of its momentum. If dM/dv and the acceleration along the 
filament are both constant, then 5 = 1/12. 

Equation (|2ip can be written in terms of scaled quan- 
tities as 

\ 10 9 Mq J \ 300 kms -1 J V 15 kpc,/ 6 V ; 

Interestingly, equation pip shows that the mechanical en- 
ergy injection rate is directly proportional to the total mass 
of the cold material. This gives a plausible theoretical expla- 
nation for the possible linear correlation betw een molecula r 
gas mass and optical luminosity in figure 9 of lEdgd (2001). 

We can use equation (|22p to estimate typical energy 
injection rates assuming that each of the parameters takes 
its typical value (e.g. Av = 300kms _1 , as seems to be the 
case in the Perseus cluster) and taking the mass of cold gas 
from column 3 of table 2. Consider the Virgo cluster; we find 
that the typical values give E ~ 5 x 10 40 ergs _1 compared 
to the observed ~ 2 x 10 40 ergs _1 . In the Perseus cluster, 
equation (I22[) gives E ~ 2 x 10 42 ergs _1 while the optical 
nebula around NGC 1275 actually emits 4.7 x 10 42 ergs _1 
in Ha and [Nil]. Thus, equation (|22p seems to provide a 
reasonable explanation for the optical luminosities of both 
systems. 

Of course, these estimates have not accounted for dif- 
ferent average filament lengths, nor the variation in relative 
velocity between the hot and cold phases in the different 
clusters. 

We can attempt to account for different average fil- 
ament lengths by comparing the sizes of th e opti cally- 
emitting regions of the clusters. iHeckman et al.l (|l989l ) give 
the sizes for Virgo, Perseus, Hydra-A, A2597 and A1795 
which have optically-emitting regions of linear diameter 14, 
>53, 13, 22 and 61 kpc, respectively. We note that the 
optically-emitting region of A1795 is dominated by a sin- 
gle filament, so that the remainder of the optical emission is 
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probably containe d within ~ 15 kp c, which provides a more 
realistic estimate. iHeckmanl (|l98ll ) also gives the size of the 
optical region in Perseus as ~ 60 kpc. Using these values we 
can estimate the average length of a filament by assuming it 
is a constant fraction of the total diameter of the optically- 
emitting region, for example 0.25. This gives average fila- 
ment lengths of ~ 4, 15, 3, 6 and 4 kpc for the filaments in 
Virgo, Perseus, Hydra-A, A2597 and A1795, respectively. 
The average filament lengths of the clusters that we do not 
have diameters for (A2199 and A478) will be taken to be 5 
kpc, which is representative of the sample. 

For the velocity shears, we have an estimate of Av ~ 
300 km s -1 along t he length of the filaments in Perseus from 
lHatch et all l|2006n , but we do not have similarly detailed ob- 
servations for any other clusters in the sample. In addition, 
linewidths at other frequencies do not give an indication of 
a single well-defined value of Av in each cluster. Therefore, 
we will assume a canonical value for Av as 300 kms -1 which 
we know applies adequately in the Perseus cluster. This will 
then be scaled to each cluster by multiplying it by the ratio 
of the cluster gravitational potential at the maximum ex- 
tent of the optically-emitting region to that of the Perseus 
cluster. That is, the value of Av that we will use for a given 
cluster will be 300 kms" 1 x (0ciustcr/<^perscus) 1//2 . El Such an 
approach is ad hoc, but does attempt to account for the dif- 
ferent environments. Henceforth we will refer to this multi- 
plicative factor by / = (^dustor/^pcrscus) 17 ' 2 , see table 3. The 
gravitational potentials are estimated from the temperature 
and density profiles for each cluster, with the assumption 
that the gas is in hydrostatic equilibriu m. Fits to the tem - 
perature and density data are taken from lPope" et alJ (|2006l ) 
and they were numerically integrated from the origin to ob- 
tain the gravitational potential at the edge of the optically 
emitting region. It was not possible to determine a gravita- 
tional potential for A2199 due to the power-law fits to the 
data. Consequently, we assigned / = 1 for this object. 

These scali ngs have b e en ap plied to the sample o f 
clusters used in iPope et all (|2006h and iPope et~ai1 l|2007h . 
The optical luminosit ies f or this sample a re ta ken from 
ICrawford et al.l l|l999h and iHeckman et all (|l989h . As can 
be seen from table 3, the expected energy injection rates are 
generally comparable with the Ha+[NII] luminosities. 

Figure 2 demonstrates a relatively good match between 
the mechanical energy injection rate, and the combined 
[NII]+ Ha optical luminosity. In fact, the best-fit shows that 
/ x 400 kms -1 provides a better description of the veloci- 
tie s. Thi s , part icularly the scaling, agrees with the findings 
of lEdgei |200ll ). It therefore seems possible that the dissi- 
pation of mechanical energy may be evident in the [Nil] 
line emission, which traces the heating rate, while there 
may be an ionising source which produces the Ha emis- 
sion. Regardless of this, the general conclusion that optical 
emission increases w ith mechanical power is consistent with 
iFerland et alJ (120081 ) who found that the H2 to H recombi- 
nation line ratios are compatible with dissipative heating. 

It is worthwhile pointing out that our estimate for 
A1795 falls below the observed optical luminosity. This is 
probably because the assumed velocity is a conservative es- 
timate given the current likely state of the system. The AGN 



1 The square-root occurs because v ~ c/) 1 / 2 . 
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Figure 2. Comparison between the mechanical energy injection 
rate and Ha + [Nil] luminosity. The diamonds show the objects 
for which the E is an upper limit. As a result we may expect 
these points to move to the left if the mass of cold gas is less than 
the upper limit. The stars denote the remainder of the sample. 
Ideally we would compare the mechanical energy injection rate to 
L([NII]), but this is not possible since we do not have L([NII]) val- 
ues for the objects in the Hcckman et al. (1989) sample. The line 
indicates where the optical luminosity is equal to the mechanical 
energy injection rate. Note that alterations in the normalisation 
could be achieved by employing a different relative velocity be- 
tween the phases. Uncertainties are shown for A478, by Crawford 
et al. (1999). The other measurements are taken from Heckman 
et al. (1989) who did not provide uncertainties. 

has been recently active and indeed has injected more en- 
ergy than any of the others in the sample, save for the one in 
Hydra-A, which is a well-known powerful AGN, so a fast out- 
flow would be expected in the central regions. In addition, 
this system would only requires An ~ 300 kms -1 to account 
for the optical emission, so it seems more than likely that 
mechanical energy can account for the optical emission. A 
similar explanation also accounts for the apparent deficit in 
the Perseus cluster. 

The question of why some clusters are brighter than 
others in their optical emission remains. The mechanical 
energy injection rate seems to be proportional to the total 
mass of cold gas so that clusters with more cold gas are likely 
to exhibit more optical emission. This leads to the question: 
why is there more cold gas in some clusters than others? The 
mass of cold gas is related to the mass flow rate, which is, in 
turn, related to the balance between heating and cooling in 
the cluster. There may also be an additional component of 
material thrown up from the central galaxy. It is important 
to note the strong effect of the velocity difference between 
the phases. This is again likely to be related to AGN activ- 
ity, so that the optical luminosity in clusters with AGN is 
likely to be larger than expected. 

4.1 Infra-red emission 

There is also significant infra- red emission from the molec- 
ular gas near cluster centres. IJaffe et al.l (|2005l ) quote lu- 
minosities of ~ 10 42 ergs~ 1 for the H2(l — 0)S(1) line in 
A2597 and suggest that this comprises roughly 1-2% of the 
total molecular NIR emission due to H2. This would imply 
a total NIR luminosity of 10 44 ergs -1 which greatly exceeds 
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Table 3. Mechanical energy injection rates calculated with equation (23) using the observed cold gas masses 
(in table 1). In column 2 we list the average filament lengths derived from the size of the optically-emitting 
region. The gravitational potential / factors are given in column 3. The estimated optical luminosities are 
shown in column 4, the observed optical luminosities are in column 5, the X-ray luminosities in column 6 are 
taken from Pope et al. (2006) and references therein, and the references are shown in column 7. Reference 3 
is Heckman et al. (1989); reference 4 is Crawford et al. (1999). Using the [Nil] fluxes given by Crawford et 
al. (1999) we have calculated L(Ha + [Nil]) rather than L(Ha). The luminosities from reference 3 are the 
sum of the Ha and [Nil] luminosities. Heckman et al.(1989) did not give uncertainties. 



Object I f E/lO^crgs- 1 L(Ho + [NII])/10 40 erg s" 1 L x /10 42 ergs" 1 Ref 



Virgo 


4 


0.8 


9 


1.9 


9.8 


3 


Perseus 


15 


1 


200 


470 


670 


3 


Hydra 


3 


0.8 


<51 


15 


250 


3 


A2597 


6 


0.8 


210 


270 


430 


3 


A1795 


4 


0,1 


<8 


82 


490 


3 


A2199 


5 


1 


<90 


5.9 


150 


4 


A478 


5 


0.5 


17 


22.7 


1220 


4 



the optical luminosity for A2597. |jaffe fc Bremerl J 19971 ) also 
give luminosities of ~ 10 40 ergs -1 and ~ 10 41 ergs -1 for the 
Ha(l - 0)S(1) line in A478 and Hydra-A, respectively. This 
implies total NIR luminosities of 10 42 — 10 43 ergs -1 respec- 
tively. Such values are comparable with the optical lumi- 
nosities of the systems, although not in the case of A2597. 
However, it is still feasible that the injection of mechanical 
energy can also account for the NIR emission. To obtain a 
mechanical energy injection rate of ~ 10 44 ergs -1 in A2597 
requires Aw ~ 700 — 800 km s -1 . Such values are compara- 
tively large, but probably close to a Mach number of unity 
for the hot gas in the cluster centre and is therefore per- 
fectly permissible. Furthermore, A2597 shows evidence for 
a rec ent, moderately powerful AGN outburst (jfJirzan et al.l 
l2004h : therefore, transonic velocities close to the cluster cen- 
tre are to be expected. The same is true for Hydra-A, which 
also shows strong NIR emission and currently hosts a strong 
AGN outburst. 



5 DISCUSSION AND SUMMARY 

The effect of mass injection into a di ffuse flow in a gravi- 
tational potential was investigated by IPittard et al. I (|2004l ) 
using ID hydrodynamical simulations. Regularly spaced 
shock-like structures were obtained under some conditions, 
which could provide an alternative explanation for the 
structures generally believed to be AGN-driven shocks. 
Other simulations revealed a differential luminosity func- 
tion which resembled those deduced from observations (e.g. 
IPeterson et~aT1l2003l ). However, along with these successes 
are a number of areas for concern. Firstly, the amount of 
material that was added was not related to the mass de- 
position rate (expected from cooling), and eventually dom- 
inated the total gas mass in the system. However, this is 
less of a problem if the flows are fountain-like. Secondly, 
the resulting X-ray luminosity was far too centrally peaked. 
This may be related to the form of the gravitational poten- 
tial used, and it would be worth revisting the model, to see 
what the effect of a realistic cluster potential, and realistic 
mass transfer rates, would have. It may also be necessary 
to include thermal conduction as well, to see whether these 



structures would be smoothed out by the same process that 
evaporates the clouds. 

Our analysis in the current work indicates that the rate 
of mechanical energy injection is comparable with the opti- 
cal luminosity in the sample. Since the [Nil] emission is a 
tracer of the heating rate of the nebular gas, it seems possible 
that the observed emission is indeed powered by mechanical 
energy injection. Furthermore, given that the energy dissi- 
pation rate is proportional to (u — v) 3 it is not surprising 
that optical emission occurs frequently in cool-core clusters, 
since AGN are frequently found in such environments. That 
is not to say that radiation from the AGN is responsible 
for the emission, but that AGN-induced bulk motions in the 
ICM are likely to cause significant interactions with the cold 
material in the central regions of clusters which may be ev- 
ident through the optical emission. It may also be revealing 
that optical emission in clusters with AGN is generally more 
luminous than would be predicted. It is also possible that 
a component of the H Q is due to the same process, though 
photoionisation also plays a role. As we have already stated, 
it is also possible that mechanical energy also powers the 
NIR emission. Again, prominent NIR emission seems to oc- 
cur in systems with powerful AGN. As a result, it is possible 
that the presence of cold material may automatically lead 
to physical processes which power the observed optical and 
infra-red emission. In this picture the injection of energy will 
necessarily be relatively constant and distributed, in accor- 
dance with observational constraints. 

This work suggests that there are plausible theoreti- 
cal reasons which could explain the rough proportionality 
between the total optical luminosity and the total mass of 
cold material observed bv lEdgell|200ll ). The total cold mass 
is likely to depend on the heating/cooling balance in the 
cluster centres, thus providing us important information on 
the deviation from equilibrium between these processes. 

Our numerical simulations (jPope et al.ll2008h indicate 
that, if formed by outflows, filaments can only be formed if 
there is a sufficiently fast wind (> few x 10 7 cms -1 ). With- 
out this wind, any cold, dense material stripped from a cloud 
would sink towards the cluster centre with the cloud, rather 
than forming a long tail behind the cloud. Thus, a slower 
wind may lead to amorphous optical emission, without fila- 
ments. Consequently, the morphology of the optically emit- 
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ting region is likely to be strongly dependent on the magni- 
tude of the relative motion between the ICM and the cold 
phase. 

To summarise, we have two main points: 

1) There is enough cold gas within the central 30 kpc, or 
so, of most galaxy clusters for its evaporation to contribute 
significantly to the density of the hot gas in this region. 
Therefore, to obtain a more complete understanding of the 
behaviour of the ICM this effect should be included in future 
numerical simulations. 

2) In general, it seems that the mechanical energy or 
momentum injection will not affect the thermal structure of 
the hot phase, but it may account for, at least, the [Nil] 
component of the optical emission of the filaments and the 
NIR emission. Under exceptional circumstances this process 
could also alter the flow dynamics of the ICM and account 
for some heating in the central regions of clusters. 
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